Unfolding cross-linkers as rheology regulators in F-actin networks 
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We report on the nonlinear mechanical properties of a statistically homogeneous, isotropic semi- 
flexible network cross- linked by polymers containing numerous small unfolding domains, such as the 
ubiquitous F-actin cross-linker Filamin. We show that the inclusion of such proteins has a dramatic 
effect on the large strain behavior of the network. Beyond a strain threshold, which depends on 
network density, the unfolding of protein domains leads to bulk shear softening. Past this critical 
strain, the network spontaneously organizes itself so that an appreciable fraction of the Filamin 
cross-linkers are at the threshold of domain unfolding. We discuss via a simple mean-field model the 
cause of this network organization and suggest that it may be the source of power-law relaxation 
observed in in vitro and in intracellular microrheology experiments. We present data which fully 
justifies our model for a simplified network architecture. 

PACS numbers: 87.16.Ka, 82.35.Rs, 62.20.Dc 



I. INTRODUCTION 

The cytoskeleton of eukaryotic cells can be described 
as a biopolymer gel or cross- 

linked network [J H0]. Thc 
principal constituent of the cytoskeleton is a stiff pro- 
tein aggregate, F-actin, that is cross-linked densely on 
the scale of its own thermal persistence length. In this 
sense the mesoscale structure of this biopolymer net- 
work differs substantially from that of synthetic poly- 
mer gels IH, which may be thought of as interconnected 
random walks. This realization of the profound differ- 
ences in mesoscale structure between the cytoskeleton, 
a prototypical scmiflcxiblc filament network, and flexible 
polymer gels calls into question the application of stan- 
dard rubber elasticity theory to describe the mechanics 
of semiflexible networks. 

A number of researchers have begun to explore the 
quantitative relation between the novel microstructure 
of semiflexible gels and their observed mechanical prop- 
erties 0, 0, 0]. More recently, it has been proposed 
that, due to their different architecture, these semi- 
flexible gels have macroscopic linear moduli that are 
generically more sensitive to their chemical composition 
@, E3, [O, El, [H, S3 than traditional flexible gels. Their 
response to point forces over mesoscopic distances is 
much more complex than suggested by the predictions 
of continuum elasticity [l^. In addition they exhibit a 
highly tunable (through network microstructure) nonlin- 
ear response to stress [lq]. 

Understanding the biophysical properties of the cy- 
toskeleton demands a closer examination not only of the 
material properties of generic semiflexible networks, but 
also of the chemical heterogeneity of the physiological 
cytoskeleton. Until recently, numerical and theoretical 



modeling of the relationship between the network ar- 
chitecture on the mesoscale and the long-length scale 
mechanical properties of these gels has focused exclu- 
sively on semiflexibility of the filaments [ItJ while ignor- 
ing the details of the chemical agents which cross-link 
them. Recent in vitro experiments, however, have shown 
that the nanoscale mechanics of the physiologically ubiq- 
uitous cross-linking proteins have a significant effect on 
the mechanics of cross- linked F-actin networks [lU . One 
particularly interesting class of cross-linkers from the 
standpoint of network mechanics are those that con- 
tain unfolding domains such as a-Actinin[l^, [2(| and 
Filamin [2ll. |22| . Such cross- linkers have protein domains 
along their backbone that unfold reversibly at a critical 
pulling force. It is still a matter of debate what the func- 
tion of these domains may be, but one may speculate 
that, by exposing new chemical groups in the unfolded 
domains, cross-linking agents such as Filamin may play a 
role in transducing local network strain into biochemical 
signals. 

In this paper we investigate the purely mechanical ef- 
fect of cross-linkers which have unfolding domains. Thc 
mechanical effect of domain unfolding occurs only at 
some finite applied stress, so the effect we wish to study 
is evident only in the nonlinear elastic response of the 
material observed under finite strain conditions. In this 
aspect the present work differs from much of the recent 
theoretical research on semiflexible network mechanics 
that focused on the linear response regime. 

We demonstrate two effects that the inclusion of un- 
folding cross-linkers has on the elastic properties of a 
otherwise generic semiflexible network. Firstly, since the 
cross-linkers can only sustain a finite maximum tensile 
stress before unfolding, their inclusion leads to shear soft- 
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cning of the network at large strain. Secondly, past the 
onset of shear softening the network spontaneously orga- 
nizes so that the population of cross-links at given tension 
grows exponentially or faster up to the unfolding force of 
the domains. Thus at moderate applied stresses the sys- 
tem appears to adjust its mechanical properties so as to 
achieve a strain state in which a significant fraction of 
its cross-linkers are poised at the unbinding transition of 
their internal domains. We refer to such cross-linkers as 
being in a critical state. In a thcrmalized system, the 
rapidly rising population of cross-links weighted towards 
the unfolding force would yield a fragile state of the ma- 
terial characterized by a broad, continuous distribution 
of relaxation time-scales via thermally excited subcriti- 
cal cross-link unfolding. This thermally excited unfolding 
may account for the broad, i.e. power law distribution 
of stress relaxation times that has been observed in in 
vivo experiments [23l l25l. I2H l27l . l39j and has been termed 
soft-glassy rheology (SGR) [28|]. While in our purely me- 
chanical model this fragile state, in which critical cross- 
linkers predominate, comes about as a consequence of 
applied stress, an in vivo network built with unfolding 
cross-linkers may generically evolve into this high sus- 
ceptibility state under the action of internal molecular 
motors {e.g. Myosin - not considered in our model). The 
necessity of applied stress to enter this fragile state in ab- 
sence of molecular motors appears to be consistent with 
the recent finding that pre-stress is necessary to replicate 
in vivo rheology measurements in artificially synthesized 
Actin-Filamin networks fill ]. 

The remainder of this article is arranged as follows. 
In section [TT] we present a simple mean-field theory for 
the mechanical effects of cross-linkers with unfolding do- 
mains. In scction imi wc discuss our numerical model used 
to study the strain response of an F-actin network with 
unfolding cross- links. Sections ITT1 and UlTl reiterate mate- 
rial presented in our initial article on this subject [2^ |. 
while the following sections go beyond the preliminary 
analysis presented therein. We present the results of this 
modeling in section HVl and compare them to our theory. 
Finally wc conclude in section [V] where we discuss the 
generality of our results and place them in the broader 
context of the mechanics of scmiflcxiblc gels and the mod- 
eling of the cytoskeleton in particular. 



II. THEORY 

We now consider how the incorporation of cross-linkers 
with unfolding immunoglobin domains will affect the 
equilibrium states of a random semiflexiblc filament net- 
work. The elasticity of the unfolding cross-linkers is 
highly nonlinear due to the presence of the many identi- 
cal folded protein domains (24 in the case of Filamin). It 
is observed at large tensile forces (~ lOOpN) that one of 
these domains will open, increasing the contour length of 
the molecule by about ~ 20nm and thereby relaxing the 
stress in the system. 



As a simple approximation, we model the force exten- 
sion curve of a single cross-linker as a sawtooth. Each 
branch of the sawtooth represents the entropic elasticity 
of a cross- linker with a fixed number of unfolded domains. 
For simplicity, we take the additional contour length gen- 
erated during any unfolding event I f to be a constant and 
the force-extension relation on each branch of the saw- 
tooth to be linear with spring constant kf. Thus the 
maximum force on each branch of the sawtooth, 
corresponds to the critical unfolding force of a domain. 
We also neglect the rate-dependence of this unfolding 
force that is found in noncquilibrium unfolding dy- 
namics. We refer to this simplified mechanical model of 
the physiological cross- linkers as a sawtooth cross-linker. 
Clearly any biopolymer cross-linker that has multiple un- 
folding domains that are accessible at physiologically rel- 
evant tensions will act in this manner. We expect the re- 
sults presented below to apply to a network cross-linked 
by any molecule of this class. 

We now consider the effects of such sawtooth cross- 
linkers on the zero frequency mechanics of an elastic net- 
work cross-linked by such agents. In the following we 
show that the considerations of force balance alone al- 
low us to relate the distribution of sawtooth cross-linker 
extensions modulo the sawtooth length to the, as yet, un- 
known distribution of local network compliances in the 
randomly cross-linked network. 

Imagine the equilibrium state of an individual pulled 
sawtooth cross- linker and the region of network surround- 
ing it both before and after a single unfolding event. We 
assume that for the relatively short sawtooth length If 
found in physiological networks, the response to one un- 
folding is in the linear regime of the surrounding net- 
work. We model the surrounding effective medium as a 
single harmonic spring with spring constant k sampled 
from some statistical distribution Pft-(fc) which encapsu- 
lates the differing local environments around each cross- 
link. Reflecting the network structure, the cross-linker 
is connected in series with the effective network spring. 
We first set the total strain on the two springs in series 
(by fixing their combined length) so that the cross-link 
is poised at the maximum extension of its current saw- 
tooth branch. Then wc consider an infinitesimal increase 
in the total extension which drives the unfolding of one 
more protein domains. 

Before the unfolding event the two springs are in force 
balance at the top of the previous sawtooth branch, so 
that kf£f = kx where x represents the extension of the 
medium spring. After the unfolding, the system achieves 
force balance on the next branch of the sawtooth force- 
extension curve so that the extension of the cross-linker is 
now increased by if — d while the extension of the medium 
spring is decreased to x — (If — d). Force balance requires 
that d, the distance between the current extension of the 
cross-linker and the edge of the next sawtooth, is given 
by d(k) = k£f/(k + kf). In other words, the combined 
system once equilibrated with the cross-linker at its max- 
imal force is now equilibrated with that cross-linker on 
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its next sawtooth branch at a smaller force. The strain 
in the surrounding medium has also decreased due to the 
extension of one more protein domain. 

To maintain force balance, the cross-linker cannot re- 
lax its length more than if — d. Upon further extension 
the cross-linker will only stretch more until another do- 
main unbinds. Thus in steady-state the sawtooth cross- 
linker will evenly sample all extensions (modulo if) be- 
tween if — d(k) and if. For a given value of the spring 
constant of the medium we expect that the extensions 
(modulo if) of the sawtooth cross-linkers Xf to be uni- 
formly distributed between the bounds given above. This 
distribution can be written as 

p(xf,k) = -^-e(x f -[i f -d(k)}), (i) 

where is a step-function. Integrating over the spring 
constant distribution Px(k), we write the probability of 
finding a given cross- linker length (modulo if) Xf\ 

P ^f) = / f PK(keff)dk ef f. (2) 

The step function fixes the lower limit on the fc-integral. 

There are several important observations to make con- 
cerning Eq. O First we have treated the unfolding events 
as being uncorrelated. In this mean field model we have 
ignored the effect of additional unfolding events in the 
surrounding medium. This simplification allows us to 
treat the surrounding network as a linear clastic clement 
with an undetermined spring constant k. We will test this 
assumption by comparing the predicted relation, Eq. [21 
to our numerical data. From this equation we note that 
the appearance of an exponential peak in the extension 
distribution P(xf,k) requires a nontrivial structure for 
the distribution of local spring constants Pjc(k). This 
distribution itself depends on the random connectivity 
and material distribution of the network. Below we will 
examine this distribution in some detail numerically and 
we defer speculations about its form for section IIV B[ 
where we discuss our numerical results. 

A second feature of Eq. [3] is that it can be rescaled by 
ij 1 in such a way that the dimensionless length distri- 
bution P(xf/if) is completely independent of the saw- 
tooth length if. Thus the distribution P (xf /If) modulo 
1 depends only on kf and Pft-(fc). We now turn to a dis- 
cussion of the numerical model we used to study random 
networks with unfolding cross- links. 

III. NUMERICAL MODEL 

We employ a simple simulational model of a statisti- 
cally homogeneous, isotropic network of semiflexible fil- 
aments in two-dimensions. These networks are formed 
in a manner identical to that of Head et ai.JU [HI, [II], 
i.e. by placing straight rods of a fixed length at ran- 
dom positions and orientations in a box with periodic 



boundaries. At points where two rods intersect a node 
is added to each rod and a cross-linker is added connect- 
ing the nodes. The cross-linkers have zero rest length. 
Rods are added until the average number of cross-links 
per rod reaches a target value that we use to parameter- 
ize the network density. A model network constructed 
by the procedure described above is shown in Figure [T] 
The sections of rod between nodes are modeled as linear 
springs with fixed clastic constant per unit length of rod. 
The cross-linkers exert no constraint torques at the nodes 
so that the rods are free to bend there with no energy 
cost. We write the Hamiltonian for each filament as 

where dl(s)/ds gives the strain or relative change in lo- 
cal contour length, and (i is the Young's modulus of 
the filament (essentially a spring constant normalized to 
1/ [length]). 

We found that all numerical minimization routines we 
tried failed to converge in the non-linear large shear 
regime when the energetics of filament bending were 
added to our simulational model. This is due to the 
large distribution of length scales in the random network, 
and the resulting poorly conditioned jacobian. At the 
same time, the nonlinear behavior of semiflexible net- 
works with freely rotating cross-links has been shown to 
be dominated by semiflexible filament stretching instead 
of bending [l7j . We therefore omit bending energies from 
our model. This omission allows for significant gains in 
computational efficiency and leads to negligible errors in 
the high-strain state, on which we focus. The effect of 
this approximation on our results will be discussed where 
appropriate. 

Similarly we anticipate that the results derived here are 
essentially independent of network dimensionality since 
network connectivity, not the dimensionality of the space 
in which the network is embedded, should control the 
collective mechanical properties of the system. 

Based on the justification given in Section!!]] we model 
the unfolding cross-linkers as non-linear springs with a 
sawtooth force extension law. Once again, the sawtooth 
branches have linear spring constant kf and length if. 
Though the physiological filament cross-linkers have a 
finite number of unfolding domains (24), we will let our 
simulated sawtooth force extension curve have an infinite 
number of branches. 

The network is sheared using Lees-Edwards bound- 
ary conditions (by adding a constant horizontal offset 
to filaments that crossed the top and bottom bound- 
aries) as shown in Figure [Hb). At the beginning of 
each strain step all nodes are moved affinely, then the 
node positions are relaxed through a conjugate gradient 
routine to a point of local force equilibrium. Since the 
cross-linker force extension curve is a sawtooth, there are 
many possible equilibrium states of the network. We wish 
to consider the adiabatic, history dependent states of a 
strained network, which would in principle require us to 



FIG. 1: (color online) Model network showing the F-actin filaments in blue and the sawtooth cross-linking agents in red. (a) 
unstrained, (b) 100% shear. 



use strain steps resulting in displacements smaller than 
the sawtooth length £f so that mechanical equilibrium is 
achieved in the earliest possible branch. To save compu- 
tational time we use a two step equilibration procedure, 
which finds a state close to adiabatic state, but allows 
for large strain steps. In the first equilibration step, we 
replace the sawtooth force law for all cross-linkers by the 
following force law: 



f = 



kfH. 



|x| < l t , 
W > If. 



(4) 



parameter 


symbol 


sim units 


scaled units 


Ig domain length 


if 


6.5 x 10" 4 


20 nm 


cross-link separation 


4 


6.6 x 10" 3 


0.2 um 


filament length 


£r 


0.2 


6 /im 


spring ratio 


*//(*«> 


0.06 - 66 


0.4 - 7 



TABLE I: Scaled values of simulation parameters. The scaled 
Ig domain length and cross-link separation were taken from 
the literature [3ll . |32| , and were the basis of the chosen sim- 
ulational values. The total filament length derives from the 
cross- link separation and the decision to use 30 cross- links per 
filament on average. This value is consistent with physiologi- 
cal values. The simulational value of spring ratios were chosen 
to be broader than than the typical physiological ratios, so as 
to capture all possible regimes. 



Thus, beyond the first sawtooth length, a cross-linking 
molecule exerts a constant contractile force of magni- 
tude kfif. The combined network of linear elastic rods 
and constant force cross-links is equilibrated. Finally, we 
reimpose a sawtooth force law for the cross-linkers and 
equilibrated the network a second time. Assuming all 
cross-links act independently as the network relaxes dur- 
ing this final equilibration step, the force on the cross-link 
must be less than kfif, so the cross-link will stay on the 
same sawtooth branch. Since the rest of the network was 
originally equilibrated at the critical pulling force, the 
sawtooth force law could not have reached force equilib- 
rium on any earlier sawtooth branch. In practice, col- 
lective relaxations of the network push individual cross- 
links onto different sawtooth branches in this final step. 
However, we found that for a variety of strain step sizes 
the quantitative behavior discussed in the results section 
was identical. This numerical technique reduced the re- 
quired computational time by a factor of ten or more, 
allowing us to simulate up to six realizations of random 
networks, each with an average of 1100 filaments and 
16500 crosslinks, for every set of material values in this 
study. With this technique, we were able to arrive at 
statistically meaningful results in a reasonable number 
of computer processor hours. 



IV. SIMULATION RESULTS 

The numerical values for all parameters in the sim- 
ulation, as well as their corresponding physical values, 
are presented in Table [H We present data for networks 
with a filament density such that there are on average 
30 cross- links per filament. The filaments have a unique 
length of Ir = 0.2 measured in units of the length of 
the unstrained, square simulation box. This gives an 
average cross-link separation of £ c w 6.6 x 10~ 3 in simu- 
lation units. For these values we find negligible system- 
size effects. The other length scale in the system is £f, 
the length of the sawtooth in our approximate form of 
the cross-linker force-extension relation. We found that 
for £r/£ c > 15 and £f/£ c constant, the network behav- 
ior is independent of £r. The distance between cross- 
links in physiolo gica l F-actin networks has been quoted 
as £ c = 0.2 /im |3lJ. For Filamin cross-links each un- 
folding domain adds if = 21 nm of length [jj^]. Thus 
the physiological ratio of sawtooth length to cross-link 
separation is £f/£ c = 0.1. We expect this ratio to be 
similar for all domain-unfolding cross-linkers. To explore 
the dependence of our results on this ratio and thus on 
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the particular type of sawtooth cross-linker, we present 
data for £f/£ c ~ 0.1 and £f/£ c ~ 0.02 (respectively, we 
used if = 6.5 x 10~ 4 and If = 1.3 X 10~ 4 in simulation 
units). 

There are two energy scales in the system determined 
by the extensional modulus of the filaments fi and the 
spring constant of the cross-linkers. To fix the energy 
scale in the problem we set the extensional modulus of 
the filaments to unity. The average spring constant for a 
filament segment can then be determined from the mean 
distance between cross-links, or, in other words, the net- 
work density: (kji) = 1/ (£ c ) = 150. Physically, actin is 
a worm-like polymer, so the effective spring constant of 
each free length of polymer should be k^ ~ ksT£p/£^ Q, 
where £ p is the persistence length of F-actin, which is ap- 
proximately 17 fj,m [3l| . Thus at room temperature, a 
0.1 fim segment of actin will have &r ~ 12 pN/nm, while 
a 0.2 yjm segment will have kfj ~ 0.7 pN/nm. Because 
of the variability in this effective spring constant as a 
function of £ c , we choose to simulate a range of actin to 
cross-linker spring constant ratios. The measured unfold- 
ing force for an Ig domain is around 150 pN [32j ~~ treat- 
ing the cross-linker as a linear spring between unfolding 
events and using the sawtooth length given above, we ar- 
rive at a physiological spring constant of kf ~ 8 pN/nm. 
Thus the relevant dimcnsionlcss ratio of spring constants 
to study for physiological systems may lie in the range be- 
tween 0.4 < kf I (kft) < 7. In simulation units, the range 
of cross-linker spring constant values reported here was 
between 10 1 < kf < 10 4 , which corresponds to a range 
of spring constant ratios between 0.06 < kf J (ka) < 66. 



A. Elastic moduli 

Since we have omitted the bending rigidity of the fila- 
ments, these networks necessarily have vanishing elastic 
moduli in their unstrained state [34|. It has been shown 
that in vitro F-actin Filamin networks do have a finite 
modulus at zero shear, though the differential modulus 
increases by a factor of 100 or more upon shearing [l8[. 
Our networks quickly develop non-zero moduli under fi- 
nite strain. To characterize the nonlinear mechanical 
properties of the network we measure the differential 
moduli daij/duij as a function of applied strain. Fig. [5] 
shows how the differential shear modulus evolves as a 
function of strain. As expected, the differential moduli 
all vanish for zero applied strain. Upon increasing ap- 
plied strain, they grow monotonically to some maximum 
where the network is stiffening most dramatically under 
further strain increments. At still larger values of the 
applied strain, the stress plateaus so that the differential 
moduli shown here decay to zero. In part (a) of Fig. [2] 
we show the differential shear modulus K(u xy ) for net- 
works with if I ' l c = 0.1 as computed from the additional 
shear extension of the network already under shear u xy . 
In part (b) of this figure we show the differential shear 
modulus K(u xy ) for networks with £f/£ c = 0.02. Though 



not shown, we also measure the uniaxial extension mod- 
ulus K for the case tf/i c = 0.02. From a comparison of 
the two measurements we determine the effective Poisson 
ratio v as function of applied strain [35| using 



1 3A^ - IK 

2 3A" + K ' 



(5) 



For intermediate deformations, where the network is 
strain stiffening, the scale of K is about four times that 
of K, so v 0.38. 

The peak in the differential modulus was observed to 
coincide with the appearance of "tears" in the network. 
These tears were actually clusters of highly extended 
cross-linkers. Eventually, a single cluster of highly ex- 
tended cross-linkers percolates across the network. Thus, 
the modulus peak occurs when a significant fraction of 
the cross-linkers are at their unfolding threshold. 

The formation of the fragile state characterized by 
cross-linker unfolding requires a high-stress state. Since 
F-actin is found to rupture under tensile loads of 
600pN [33], one may ask whether F-actin rupture pre- 
cludes the formation of the fragile state characterized by 
cross-linker unfolding. To check this we measured the 
tension along the actin filaments at the onset of cross- 
linker unfolding; we found that, although there are multi- 
ple cross- linkers connected to each filament, the force due 
to adjacent cross-linkers tend to cancel, so that stresses 
do not accumulate along the actin filaments. Using a 
conservatively high estimate unfolding force of Filamin, 
150 pN, we find that, at the onset of cross-linker un- 
folding, the stresses on individual segments of actin fil- 
ament (between cross-links) were all below 4 times the 
cross-linker unfolding tension (600pN), and 90% of fila- 
ments segments were at tensions below 3 times that value 
(450pN). At larger strains these forces increase, and at 
strains 1.5 times the differential modulus peak, up to 20% 
of the actin filaments may be under more than 600pN of 
tension, effectively breaking those strands. The break- 
down of the network at large strains is safely above the 
region of interest in this work making the formation of 
the fragile state possible. Finally, we note that the un- 
folding tension of the cross-linker is actually an out-of- 
cquilibrium quantity. At slower rates of tensile loading 
the cross-linkers will unfold typically at lower tensions 
due to thermally activated processes. Such thermal ef- 
fects at lower loading rates, serve to widen the range of 
strain states in which cross-linker unfolding occurs before 
F-actin failure. We will discuss the scaled tensions in the 
network further in Section fVl 

Since the tearing of the network corresponds to a seg- 
regation into actin-rich and cross-linker rich regions, we 
can thus approximate the compliance of the network as 
the compliance of a composite system. One part of the 
composite is the filament dominated parts of the network 
in which the cross-linkers are not greatly extended and 
the other part is the region of large cross-linker stretch. 
The effective modulus of the composite system can then 
be approximated as two nonlinear springs in series. 
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FIG. 2: (color online) Differential moduli versus strain for several different cross-linker spring constants kf. Modulus K is 
given in simulation units, (a) Shear modulus versus shear strain for £f/£ c = 0.1. (b) Shear modulus versus shear strain for 
£f/£ c = 0.02. (c) Shear modulus versus shear strain, for constant cross-linker force / = kf£f with £f/£ c = 0.02. From lowest 
to highest curves in each graph, the cross-linker spring constants are respectively 10, 100, 200, 400, 600, 800, 1000, 2000, 4000, 
6000, 8000, 10000 and inextensible cross-links. (From Section ITVl the average filament segment spring constant in these units is 
(/c_r) « 150.) (d) Strain value u™y X corresponding to the peak differential modulus from the data in (c) (lower points) and (a) 
(upper points) as a function of kf. The curved lines are fits to v/fe/. The straight line is a fit to ax + b for values of kf > 1000 
and £f/£ c = 0.1. 



The force law for the filament dominated parts of the 
network can be inferred from the top curve in Fig. [2Jb), 
which shows the differential shear modulus for networks 
with non-compliant cross- links. The shear modulus in- 
creases steadily at low shear as the individual filaments 
align with the shear direction. Above 200% strain, the 
modulus starts to plateau. In the limit of complete shear 
alignment, the extension of individual filaments is linear 
in the shear and thus the differential modulus is constant, 
with value K (u xy ) = n ave ku, where n ave is the average 
number of filaments per unit cross sectional length. At 
low and intermediate values of the shear, the dependence 
of the modulus on absolute shear is complicated, reflect- 
ing simultaneous filament extension and alignment. In 
practice, we observe the growth of the modulus in this 
regime to be approximately linear in the shear. 

In contrast, the part of the network dominated by 
highly extended cross-links can only sustain a finite maxi- 
mum stress per unit area, since each cross-linker can only 



exert a maximum force oikftf. In the limit of complete 
shear alignment, the maximum stress is given by 

O'max = n cr kf£f, (6) 

where f~L f f*f IS the number of cross-linkers per unit of cross- 
sectional length that take part in the crazing, or splitting, 
of the network. f~L f f*f IS determined by the shortest percola- 
tion path across the network. We surmise that the differ- 
ential modulus of the composite system reaches it maxi- 
mum value at the onset of crazing. The expression above 
for the maximum stress implies that the onset of crazing 
is determined purely by the product kftf. Indeed, Fig- 
ure [He) shows the modulus for sheared networks where 
the cross- linker force law is simply that given in Eq.[5] es- 
sentially a constant force cross-link with force kflf. The 
curve is nearly identical to Figure [2fb), for the sawtooth 
force law. A constant force law is by definition indepen- 
dent of network extension, so the differential modulus 
ddij/duij is zero in a stretched, cross- linker dominated 
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region of the network which has plateaued at the maxi- 
mum StreSS (Jmax- 

We now consider shear strain that maximizes the dif- 
ferential shear modulus. The maximum differential shear 



modulus occurs at the crossover point. 



"•xy 1 

applied stress approaches n cr kf£f as shown in Eq. [6] 
For smaller kf/kn, the crossover occurs in the linear 
growth region of the differential modulus, where o~ xy ~ 
n a vekRu1y, so the shear at maximum scales as u™y X ~ 
y/kf/kR. For larger kf/kji the crossover occurs in the re- 
gion of constant differential shear modulus, so the shear 

kf/kft. The plots in Fig- 
versus kf for the data in Figure [2a) 
and (b). The data for lfji c — 0.02 is consistent with 
kf. For if I in = O.f, the strain at maximum 
■ kf < 1000, then approaches a 
linear dependence on kf at higher kf. 

For larger relative values of if the splitting of the net- 
work, and therefore the drop-off of the differential mod- 
uli, is suppressed at lower shear values. It does not oc- 
cur until a finite fraction of the cross-linker population 
is stretched beyond its initial sawtooth length if. Still, 
in the limit of large strain, the modulus is mainly deter- 
mined by the combination F max = kfif. In the next sec- 
tion we will show that most stretched cross-linkers reach 
equilibrium near F = F max , so it is natural that the net- 
work elastic response is essentially that of a network with 
constant force cross-linkers. 

In living cells, the action of molecular motors might 
lower the threshold of cross-link extension beyond the 
sawtooth length by prestressing the network. Indeed, 
it has been shown [181 ] that prestress is necessary to 
replicate physiological conditions in in vitro F-actin Fil- 
amin systems. Physiologically, unfolding-based behavior 
should manifest itself once the average stress per filament 
exceeds the Ig domain unfolding stress of ps 150 pN [32| . 
For a dense network with 10 filaments per jim 2 of cross 
section the critical prestress would amount to 1.5 kPa, 
while for sparse network with on the order of 1 filament 
per /im 2 the critical prestress would be 150 Pa or less. 



when the 



at maximum scales as u™y X 
urc[2jd) shows u mc 



"xy - V ' /' 

scales as u ml 



B. Cross-linker extension statistics 



For values of kf/{kn) < 10 the statistical weight for 
finding a cross-link extension (modulo if) appears expo- 
nentially enhanced towards length i f where the domains 
unbind. For values kf/(kn) > 10 the statistical weight 
for finding a cross- link extension (modulo if) grows faster 
than exponentially near length if. For kf / (k^) > 4 the 
distributions of cross-link lengths were identical for both 
measured values of if . At lower k f there were significant 
differences between the if/£ c — 0.1 and £f/i c = 0.02 
data, with the former showing a less pronounced pilcup 
towards length if. The source of this deviation could 
be non-linearities in the response of the softer networks 
over the longer length scale of if/i c = 0.1 ~ our effective 
spring model assumed linear response. 

According to the model developed in Section [Til the 
cross-link length distribution is determined completely 
by the distribution of local effective network spring con- 
stants. In particular, there is nothing in Eq.[^]that would 
create an exponential in the length distribution unless 
that exponential were already contained in Px(k), the 
effective spring constant distribution function. 

Figure (2Kb) shows the measured distribution of ef- 
fective network spring constants experienced by the 
cross-linkers sampled in Figure [3ja). This distribu- 
tion was found by probing one cross-link at a time, 
changing the force on the cross-link and numerically re- 
minimizing over lattice displacements to find the corre- 
sponding change in length. We find that there is indeed a 
non-trivial and rapidly decaying distribution of effective 
spring constants. For high kf the spring constant dis- 
tribution approaches a power-law with exponent —2.25, 
while for lower kf the distribution is truncated at high 
effective k, becoming nearly exponential. 

We may use Eq.[5]to calculate P(xf), the expected dis- 
tribution of cross- linker lengths modulo I f , for the cases 
of exponential or power law Px(k). For the the power 
law case with Pk^) ~ k~ p , Eq. [5] yields a cross-link 
length distribution of the form 



L.1-P 

K f 



p(p-l)£ f 



(P-1) 



p 



f-i - x f 
■'■/ 



l f-Xf 



p-1' 



(7) 



The crazing or splitting of the network discussed in 
the last section occurs once the local strain forces expe- 
rienced by cross- links exceeds kfif, the maximum force 
sustainable by a cross-link sawtooth branch. We found 
that, in the subset of all cross-links which had been ex- 
tended through at least one "unfolding event" (i.e., those 
cross-links which had changed sawtooth branches at least 
once) , a characteristic distribution of equilibrium lengths 
(modulo if) emerges that is independent of total strain. 
Figure[3]Ja) shows the measured equilibrium distributions 
of cross-link lengths, modulo the sawtooth length if, for 
a representative set of strained networks with if ji c = 0.1 
or if ji c = 0.02 and several values of spring constant kf. 



This form for P(xf) diverges at if for any p > 0. Alter- 
nately, if the effective spring constants were exponentially 
distributed, with form Px(k) ~ e~ k / k , then P(xf) would 
take the form 



■ cxp 



k f( x f 



kXf 



where the constant k is an undetermined material pa- 
rameter and r is the incomplete Gamma function. This 
expression grows exponentially near Xf =lf. 
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FIG. 3: (color online) (a) Distribution of normalized cross-linker lengths £/£f modulo 1 in equilibrated networks with, from 
shallowest to steepest slopes respectively, kf = 100, 600, 2000, and 6000 (where (&r) = 150). For the three smaller values of 
kf, values were given for £f/£ c = 0.1 (crosses) and for £f/£ c = 0.02 (boxes). For kf = 6000 no significant pileup was measured 
for £f/£ c — 0.1. In plots (b) - (d) we only plot values for £f/£ c — 0.02. (b) Measured effective medium spring constant 
distributions for kf = 600, 1000, 2000 and 6000 (higher kf have longer tails). The dashed line shows fc -2 ' 25 . (c) Distribution of 
normalized cross-linker lengths £/£f modulo 1 for kf = 600 (boxes) compared to distribution predicted by Eq. [3]for the data in 
(b) (crosses), (d) Distribution of normalized cross-linker lengths £/£f modulo 1 for kf = 6000 (boxes) compared to distribution 
predicted by Eq. [2] for the data in (b) (crosses). 



Eqs. [7] and [5J together with the data in Figure EJb), 
are consistent with our observation that the cross-linker 
length distribution grows exponentially for low kf and 
even faster for high kf, since the respective forms for 
Pif(fc) for these two cases are exponential and power law. 
Figure GUc) shows the result of numerically evaluating 
Eq.[2]on the data from Figure[3Kb) for kf / (k R ) = 4. The 
agreement between the theoretical prediction for P(xf) 
and the actual measured distribution is very good at high 
Xf. Similarly, Figure [3Jd) shows the result of numerically 
evaluating Eq. [5] for the measured Pft-(fc) for kf/ (ku) = 
40. The agreement between theoretical prediction and 
measurement is outstanding. 

Figure g] again shows P/f(fc) for kf/{k R ) = 40. The 
distribution function for sawtooth force law cross-linkers 
is compared to the distribution for a network which 
was created and equilibrated under strain with linear 
force law cross-linkers. Both distributions are identical, 
demonstrating that the presence of the sawtooth force 



law does not determine the distribution of Pr-(/c). This 
justifies our use of a mean-field approach in Section HT1 
since Px(k) is not influenced by P{xf). The further 
implication is that the origin of the measured Pr-(&) is 
purely geometrical. 

We have seen that for large values of the sawtooth 
spring constant kf the local spring constant distribution 
in the random network has a power law decay in the 
large k regime as shown by the boxed points in Fig. \3[b) . 
These points are fit by the dashed line representing a 
power law with exponent —9/4. The power law behavior 
is also clearly visible in the log-log plot in Fig. [4] which 
shows the higher kf data with £f/£ c = 0.02. For smaller 
values of kf, e.g. the crosses in the Fig. [3Jb), the local 
spring constant distribution exhibits a clearly exponen- 
tial decay. Considering the full set of the numerical data 
shown in Fig. [3Jb) for varying values of kf ranging (in 
simulation units) from 600 to 6000, we note the follow- 
ing trends. First, the distribution appears to transition 
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FIG. 4: (color online) Measured effective medium spring con- 
stant distributions for kf/ (kn) = 40. Box points are for a 
network with sawtooth cross-linkers, while the cross points 
are for a network with linear spring cross-links. Plots for 
lf/l c = 0.02. 

from a power law form at higher k in the case of stiffcr 
cross-linker springs, i.e. larger values of kf, to an expo- 
nential form for softer cross-linker springs. Secondly, the 
maximum values of local spring constants observed in the 
system appear to be generically smaller than the cross- 
linker spring constant for that system. Finally, each data 
set for a given kf appears to follow the power law decay 
for smaller values of k before reaching an exponentially 
decaying terminal regime at the largest measured values 
of k. 

We now speculate as to the origin of these phenomena 
beginning with an examination of the power lay decay 
seen for smaller values of k in all the data sets. Since 
this effect is more robust at larger values of the cross- 
linker spring constants, consider first the limiting case of 
perfectly incompliant cross-linkers. Here the local effec- 
tive spring constant of the network must be completely 
determined by the filaments, and, in particular, by the 
distribution of filament lengths between cross- links. The 
distances between cross-links, £ c , for a random network 
such as ours are exponentially distributed, with 

P(l c )d£ c ~ e-^«>d£ c . (9) 

as shown by Kallmes and Corte [36[. Since the spring 
constant of each filament segment is proportional to l/£ c , 
the probability of the occurrence of a filament segment 
spring constant between k and k + dk takes the form 

P(k)dk ~ e- 1 /«^>fe) fc -2 dA ._ ( 10 ) 

Furthermore, if we assume that softer regions of the 
network from which the smaller values of local spring 
constant k are measured form small pockets in a mate- 
rial that is generally less compliant, then it is reasonable 
to suppose that the deformation imposed on the network 
due to the external forces applied in those softer regions 
will lead to a highly localized deformation involving only 



a few network springs. Then the local effective spring 
constant around those cross-links is determined primar- 
ily by a few filament segments whose effective spring con- 
stants statistics will be governed by Eq. |TD] For large 
values of k, this expression approaches k~ 2 , which is con- 
sistent with our findings for large kf . The exponential in 
Eq. [TU] becomes unimportant for 1 < (i c ) k so that for 
150 < k we expect to see a only power law distribution 
of local spring constants with an exponent of —2. 

The effective spring constant of any two springs in se- 
ries is less than the spring constant of either of them. 
Since all filaments are connected to one another through 
cross-linkers with spring constant kf, it is clear that this 
power law tail in Eq [TUJ cannot extend to values of k 
greater than kf. If we consider the large- k tail of the 
spring constant distribution, we must look at rarely found 
regions in the network that lie on chains of anomalously 
stiff network springs corresponding to very short network 
filaments. Each of these chains of springs is made up of 
a number of statistically independent springs connected 
in series. In order to find an extremely large value of 
the effective spring constant k it must be that for one of 
the force paths all of the constituent spring constants are 
large, since the compliance of the springs in series will 
be dominated by any single soft spring. We expect the 
probability of such a rare event to be Poisson distributed 
so that, in the high-A; tail, the distribution Pxik) takes 
the form 

P K (k) ~ H{k)e~ k / k (11) 

where H is some regular function characterizing the 
small-A; behavior of the distribution (H(x) — » const as 
x — > oo) discussed above and the constant k is undeter- 
mined by this heuristic argument. Such an exponential 
distribution is consistent with the data for low kf pre- 
sented in Figure ^[b). 

Interestingly, we do indeed see that each distribution 
has an exponential tail for values of k approaching an 
upper limiting value of kf . The value of this upper limit 
appears to be oc kf with a coefficient of proportionality 
of order unity. This suggests that if the stiff chains of 
springs associated with the stiffest region of the material 
contains of order n springs in series then there must be of 
order n such chains in parallel to account for the upper 
limit. 

The principal effect of changing kf on the elasticity 
of the random network is to push the cross-over from a 
power law decay to the exponential tail to larger values 
of k. We have presented arguments for both this lower- A: 
power law and higher-fc exponential dependence of the 
local spring constant distribution. We cannot account 
for approximately ten percent discrepancy between the 
measured and predicted exponents for the power law and 
we cannot estimate using these arguments the cross-over 
value of k between these two behaviors. Nevertheless we 
see regardless of the form of the local spring constant dis- 
tribution, Eq.[5]accurately relates that distribution to the 
cross-linker extension distribution as shown by Fig. [21c) 
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and Fig. [3td) 



DISCUSSION 



In this work we have studied a composite network 
of linear elastic elements having randomly distributed 
spring constants cross-linked by linkers that have the 
sawtooth force-extension relation common to proteins 
with repeated unfolding domains. While this system is 
clearly an oversimplification of both the chemical com- 
plexity and semificxiblc character of the cytoskeleton, 
these networks retain one important micro-architectural 
feature of the F-actin networks in that forces must 
propagate from filament to filament through a linking 
molecule exhibiting a highly nonlinear (sawtooth) force 
response to strain. We suspect that our model system 
can thus inform the understanding of the mechanical 
properties of physiological cytoskeletal networks, which 
are the subject of recent theoretical work and are be- 
ing probed experimentally with increasing quantitative 
accuracy [23J, |37|, |38| . 

Our most striking result is the observation of the de- 
velopment of a highly fragile mechanical state in the net- 
work at large strains. At moderate strain, a nontrivial 
number of cross-linking molecules reach a critical state 
where they are poised to unfold another domain. The 
presence of fluctuating internal stresses in the cytoskele- 
ton produced by variations in molecular motor activity 
and/or thermally generated fluctuations can act on this 
highly fragile state to produce strain fluctuations at all 
frequency scales, due to the broad distribution of local 
energy wells in the system. Thus, the formation of this 
critical state under applied stress may explain a partic- 
ular feature of the low-frequency strain fluctuations as 
observed by intracellular microrheology. 

The inclusion of sawtooth cross-linkers in the analy- 
sis of the mechanics of semificxiblc networks may ex- 
plain some features of the nonlinear elasticity of the cy- 
toskeleton. It is well-known that in vitro F-actin net- 
works when cross-linked by rigid, inextensible molecules 
generically show strain hardening [l6| at least at network 



densities consistent with the affine response regime [12 



We did not include the nonlinearity of the longitudinal 
compliance of the actin filaments, which accounts for 
this strain-hardening. The mechanics of unfolding cross- 
linkers presents an additional nonlinear effect that de- 
creases the differential modulus. Taking both effects into 
account, we expect that highly strained F-actin/Filamin 
gels to strain harden less than than that predicted by 
RefQ] and eventually show a strain softening regime at 
high strains. 

While strain hardening is a hallmark of highly cross- 
linked semiflcxible networks^the in vivo cytoskeleton has 
been seen to strain soften[39(. Based on this work we 
suggest that strain softening in densely cross- linked semi- 
flexible networks may be the result of the unfolding of 
specialized cross- linking proteins in it. From our numer- 



ical data it is clear that cross-linking of these networks 
by unfolding linkers can generate a strain softening ma- 
terial at large enough stresses so that domain unfolding 
occurs. It remains to be seen how or whether the interac- 
tion of these cross-linkers with motor proteins shifts this 
nonlinear elastic regime to lower applied strains. 

We believe that this work suggests the appropriate 
theoretical framework for understanding the underlying 
mechanism by which the system reaches this highly frag- 
ile state. We note, however, that much remains to be 
done in order to develop this understanding into a the- 
ory that makes quantitatively accurate predictions. In 
addition it is clear that more precise numerical explo- 
rations of the network are required in order to better 
characterize both the local elastic constant distribution 
in the network as well as the cross-linker length distri- 
bution. Nevertheless, it appears that the creation of this 
fragile state is a robust phenomenon. 

In addition, recent experiments by the Janmey group 
report that rigidly cross linked semiflexible networks 
generically show a large, negative first normal stress 
coefficient pol] . We know of no such measurements on the 
in vivo cytoskeleton. Noting that sawtooth cross-linking 
of these networks changes their nonlinear elasticity, we 
believe it will be important to study the effect of un- 
folding cross-linkers on this nonlinear elastic response as 
well. 

There are a number of other extensions of this work 
that remain to be considered. We are currently work- 
ing to add thermalized subcritical unfolding, along with 
energy input through the action of simulated molecular 
motors. The non-equilibrium steady state thus estab- 
lished will be closer to the conditions inside cells. The 
strain rate dependent unfolding force (30j then depends 
on the frequency spectrum of molecular motor activity, so 
that the effective unfolding threshold, and thus all other 
forces in the network, may be shifted down by a factor of 
2 or more. The development of a more complete model 
that includes the effect of internally generated random 
stresses due to the action of molecular motors will be 
an important step towards the direct calculation of the 
low-frequency dynamics of this biopolymcr gel of funda- 
mental biological importance. 



Acknowledgments 

BD and AJL thank J.C. Crocker for providing unpub- 
lished data and for enjoyable discussions. BD also thanks 
David Morse for enlightening discussions. AJL was sup- 
ported in part by nsf-dmr0354113 BD acknowledges the 
hospitality of the UCLA department of Chemistry and 
the California Nanoscience Institute where part of this 
work was done. BD also acknowledges partial support 
from[nsf-dmr0354113 and the Institute for Mathematics 
and its Applications with funds provided by the National 
Science Foundation. 



11 



[1] B. Alberts, D. Bray, J. Lewis, M. Raff, K. Roberts, and 
J.D. Watson, Molecular Biology of the Cell, 3rd edition 
(Garland, New York, 1994); P. A. Janmey Curr. Opin. 
Cell Biol. 3, 4 (1991). 

[2] T.D. Pollard and J. A. Cooper, Ann. Rev. Biochem. 55 
987 (1986). 

[3] E.L. Elson Annu. Rev. Biophys. Biochem. 17, 397 (1988). 
[4] M. Rubenstein and R.H. Colby, Polymer Physics (Oxford 

University Press, London, 2003). 
[5] P.A. Janmey, S. Hvidt, J. Lamb, T.P. Stossel Nature 345, 

89 (1990). 

[6] K. Kroy and E. Frey, Phys. Rev. Lett. 77, 306 (1996). 
[7] R. Satcher and C. Dewey, Biophys. J. 71 109 (1996). 
[8] F.C. MacKintosh, J. Kas and P.A. Janmey, Phys. Rev. 

Lett. 75, 4425 (1995). 
[9] D.A. Head, A.J. Levine and F.C. MacKintosh, Phys. 

Rev. Lett 91, 108102 (2003). 
[10] J. Wilhelm and E. Frey, Phys. Rev. Lett. 91, 108103 

(2003). 

[11] D.A. Head, F.C. MacKintosh and A.J. Levine, Phys. 

Rev. E 68, 025101 (R) (2003). 
[12] D.A. Head, F.C. MacKintosh and A.J. Levine, Phys. 

Rev. E 68, 061907 (2003). 
[13] Alex J. Levine, D.A. Head, and F.C. MacKintosh J. 

Phys.: Condens. Mat. 16, S2079 (2004). 
[14] B.A. DiDonna and T.C. Lubensky Phys. Rev. E 72 

066619 (2005). 

[15] D.A. Head, A.J. Levine, and F.C. MacKintosh Phys. 
Rev. E Accepted for publication (2005). 

[16] M.L. Gardel, J.H. Shin, F.C. MacKintosh, L. Mahade- 
van, P. Matsudaira, D.A. Weitz Science 304, 1301 (2004). 

[17] P.R. Onck, T. Koeman, T. van Dillen, and E. van der 
Giessen Phys. Rev. Lett. 95 178102 (2005). 

[18] M. L. Gardel, F. Nakamura, J. H. Hartwig, J. C. Crocker, 
T. P. Stossel, and D. A. Weitz, Proceedings of the Na- 
tional Academy of Sciences 103 1762-1767 (2006). 

[19] S. Labeit and B. Kolmerer Science 270, 236 (1995). 

[20] M. Rief, M. Gautel, F. Oesterhelt, J.M. Fernandez, H.E. 
Gaub Science 276, 1109 (1997). 

[21] I. Schwaiger, A. Kardinal, M. Schleicher, A. A. Noegel, 
and M. Reif Nature Struct. Mol. Biol. 11, 81 (2004). 

[22] D.J. Brockwell, G.S. Beddard, E. Paci, D.K. West, PD. 



Olmsted D.A. Smith, and S.E. Radford Biophys. J. 89, 
506 (2005). 

[23] B.D. Hoffman, G. Massiera, and J.C. Crocker 

|cond-mat/0504051| (2005) . 
[39] M.F. Coughlin, M. Puig-de-Morales, P. Bursac, M. 

Mellema, E. Millet, and J.J. Fredberg Biophys. J. 90, 

2199 (2006). 

[25] B. Fabry, G.N. Maksym, J.P. Butler, M. Glogauer, D. 
Navajas, and J.J. Fredberg, Phys. Rev. Lett. 87, 148102 
(2001). 

[26] R.E. Laudadio, E.J. Millet, B. Fabry, S.S. An, J.P. But- 
ler, and J.J. Fredberg, Am J Physiol Cell Physiol 289, 
C1388 (2005). 

[27] P. Bursac, G. Lenormand, B. Fabry, M. Oliver, D.A. 

Weitz, V. Viasnoff, J.P. Butler, J.J. Fredberg, Nat. Mat. 

4, 557 (2005). 
[28] P. Sollich, Phys Rev E, 58, 738 (1998). 
[29] B.A. DiDonna and Alex J. Levine, Phys. Rev. Lett. 97 

068104 (2006). 

[30] E. Evans and K. Ritchie Biophys. J. 72, 1541 (1997). 
[31] E. Frey, K. Kroy and J. Wilhelm, in Dynamical Networks 

in Physics and Biology,ed. D. Beysens and G. Forgacs 

(EDP Sciences Springer, Berlin, 1998). 
[32] Furuike, S., T. Ito, and M. Yamazaki, Febs Letters 498, 

72 (2001). 

[33] Y. Tsuda, H. Yasutake, A. Ishijima, and T. Yanagida, 
Proceedings of the National Academy of Sciences 93 
12937, 1996 

[34] M. Kellomaki, J. Astrom, and J. Timonen Phys. Rev. 
Lett. 77, 2730 (1996). 

[35] L.D. Landau and E.M. Lifshitz Theory of Elasticity 3 rd 
Ed. Pergamon Press, Oxford (1986). 

[36] O. Kallmes and H. Corte, Tappi 43, 737 (1960). 

[37] A.W.C. Lau, B.D. Hoffman, A. Davies, J.C. Crocker, and 
T.C. Lubensky Phys. Rev. Lett. 91, 198101 (2003). 

[38] L. Deng, X. Trepat, J.P. Butler, E. Millet, K.G. Morgan, 
D.A. Weitz, and J.J. Fredberg Nat. Mat. 5, 636 (2006). 

[39] J.J. Fredberg Private Communication (2006). 

[40] F.C. MacKintosh and P.A. Janmey Private Communica- 
tion (2006). 



